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Abstract. We investigate theoretically the low-temperature physics of a two- 
component ultracold mixture of bosons and fermions in disordered optical lattices. 
We focus on the strongly correlated regime. We show that, under specific conditions, 
composite fermions, made of one fermion plus one bosonic hole, form. The composite 
picture is used to derive an effective Hamiltonian whose parameters can be controlled 
via the boson-boson and the boson-fermion interactions, the tunneling terms and the 
inhomogeneities. We finally investigate the quantum phase diagram of the composite 
fermions and we show that it corresponds to the formation of Fermi glasses, spin 
glasses, and quantum percolation regimes. 
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1. Introduction to disordered quantum systems 

1.1. From condensed matter physics to ultracold atomic gases 

Quantum disordered systems is a very active researcli field in condensed matter physics 
(CM) initiated by the work by P.W. Anderson P who first pointed out that quenched 
[i.e. time-independent) disorder can dramatically change the properties of a quantum 
system compared to its long-range ordered counterpart. Hence, disorder plays a 
central role in modern solid state physics as it can significantly alter electronic normal 
conductivity 2J, superconductivity and the magnetic [HE] properties of dirty alloys. 

It is by now clear that studying disordered systems is comparable to opening the 
Pandora box. On the one hand, disorder introduces a list of non-negligible difficulties. 
First, averaging over disorder usually turns out to be a complex task that requires 
either original methods such as the replica trick and supersymmetry [H] or numerical 
computations using huge samples or a large number of repetitions with different 
configurations. Second, the possible existence of a huge number of excited states with 
infinitely small excitation energies leads to complex quantum phases such as bose glasses 
[Zj or spin glasses Third, the interplay of kinetic energy, particle-particle interactions 
and disorder is usually a non-trivial problem [S] which is still challenging. On the other 
hand, disorder leads to an extraordinary variety of physical phenomena such as, for 
example, Anderson localization P, quantum percolation jH], or quantum frustration 

m- 

Studies of quantum disorder in CM have several limitations. First, the disorder 
cannot be controlled as it is fixed by the specific realization of the sample. In particular, 
one cannot switch adiabatically from one configuration to another. Second, the particles 
are fermions (electrons) and the inter-particle interaction corresponds to the Coulomb 
long-range potential. Third, theoretical studies rely on toy models and experiments do 
not provide control parameters. Thus, it would be highly desirable to consider new kinds 
of disordered quantum systems. As shown in this work and also discussed in other papers 
[m 1121 , one exciting possibility is that of ultracold atomic gases where these problems 
can be revisited with unique parameters control and measurements possibilities. 

1.2. Ultracold atomic gases 

Following the recent progress in cooling and trapping of neutral atoms ^3], dilute atomic 
Bose- Einstein condensates (BEC) jT3], degenerate Fermi gases (DFG) ^21 and mixtures 
of both jin] are now currently produced at the laboratory. In these systems almost all 
parameters prove to be highly controllable J7j: (i) Using standard techniques, one can 
easily control the trapping potential, the size, the density and the temperature of the 
atomic gas. (ii) Due to the strong dilution, the contact interactions are usually small 
and can be computed ab initio. In addition, their sign and strength can be controlled 
using Feshbach resonances |EJ- (iii) Quantum statistics is also a degree of freedom 
as one can use bosons or fermions. (iv) Finally, inhomogeneous trapping potentials 
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can be designed almost at will using standard optical techniques. In particular, periodic 
potentials (optical lattices) with no defect nor phonons can be designed in a wide variety 
of geometries [TH] . Controlled disordered [20] or quasi-disordered [211 122] potentials can 
also be optically produced opening new possibilities. Hence, equilibrium and transport 
properties of interacting Bose-Einstein condensates in random potentials have been 
investigated in recent experiments |2S1I211I2S1- 

In recent works ^21 121]; we have shown that Fermi-Bose mixtures in an optical 
lattice with diagonal {i.e. with random on-site energy) randomness constitute 
a case study of quantum disordered systems. In the strongly correlated regime 
(strong interactions), the dynamics reduces to a universal effective Hamiltonian whose 
parameters can be controlled through the boson-boson and fermion-boson interactions, 
or the depth of the periodic and disordered potentials. It results in a rich quantum 
phase diagram ranging from Fermi glasses and Fermi Anderson localization to quantum 
spin glasses and quantum percolation [I^. 

This paper reviews our works on strongly correlated Fermi-Bose mixtures in 
disordered optical lattices. It further details the results discussed in Ref. [12]. In 
Ref. 126], we have presented a complete study of the system and discussed a wide 
variety of cases, including several composite types and quantum phases. As in such a 
rich system, completeness competes with conciseness, we think that a shorter review 
would be useful for non-specialist readers. It is the aim of this paper to provide such a 
concise review of our findings. We have thus chosen to restrict the present analysis to 
a specific case that proves to be non the less paradigmatic but also one of the richest 
ones. For more details, the interested reader should refer to Ref. [2E]- 

The paper is organized as follows. In section [21 we introduce the model describing 
the considered system (12.11) and the composite fermions formalism that leads to the 
effective hamiltonian ()2.2I and 12. 3^ . Section [S] shows the results for the non-disordered 
and weakly disordered lattices and in section |3] the strong disorder limit, the spin glass 
limit, is discussed. Finally, we summarize and discuss our results in section |H1 

2. Strongly correlated composite fermions in inhomogeneous optical lattices 
2.1. The Fermi-Bose Hubbard Hamiltonian 

Consider a mixture of ultracold bosons (b) and spinless (or spin-polarized) fermions 
(f) trapped in an optical lattice. In addition, the mixture is subjected to on-site 
inhomogeneities consisting in a harmonic confining potential and/or in diagonal disorder. 
Eventually, the periodic potential and the inhomogeneities are different for the bosons 
and the fermions. However, we assume that the lattices for the two species have the 
same periodicity. The lattice sites are indexed by i and are thus the same for the 
fermions and the bosons. In all cases considered below, the temperature is assumed to 
be low enough and the potential wells deep enough so that only quantum states in the 
fundamental Bloch bands for both the bosons and the fermions are populated. Notice 
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that, this requires that the filhng factor for fermions is smaller than 1 i.e. the total 
number of fermions is smaller than the total number of lattice sites A^. 

We use the Wannier basis of the fundamental Bloch band which corresponds to 
wave-functions well localized in each lattice site j2Zj for both the bosons and the 
fermions. This basis is particularly well suited for the strongly correlated regime that 
is investigated here |^. Then, the second quantization Hamiltonian reduces to the 
Fermi-Bose Hubbard (FBH) model El 122] : 

%BH = - E [^b^l^. + ^f///.- + ^-c- 

(ij) 

'V 



i 

+E 



^■Ui^Ui - 1) + Uuinii 



1 r 

-Hi Ui - HiTni 



(1) 



where bi and are bosonic and fermionic annihilation operators of a particle in the 
i-th site and Ui = blbi, rui = fjfi are the corresponding on-site number operators. The 
FBH model describes: (i) nearest neighbor boson (fermion) hopping, with an associated 
negative energy, — Jj^ (~^f)i (ii) on-site boson-boson interactions with an energy V, 
which is supposed to be positive {i.e. repulsive) in the reminder of the paper; (iii) 
on-site boson-fermion interactions with an energy U, which is positive (negative) for 
repulsive (attractive) interactions; (iv) on-site energy due to interactions with a possibly 
inhomogeneous potential, with energies —/if and —h]. Eventually, — /xf and —hj also 
contain the chemical potentials in grand canonical description. 

In the following, we investigate the properties of the strong coupling regime, i. e. 
V,U ^ J]g f and we derive a one-species effective Hamiltonian using a perturbative 
development up to second order in J]^ f/V. 

2.2. Zeroth- order perturbation: formation of fermion composites 

In the limit of a vanishing hopping (Jj^ = Jf = 0) with finite repulsive boson-boson 
interaction V, and in the absence of interactions between bosons and fermions {U = 0), 
the bosons are pinned in the lattice sites with exactly hi = + 1 bosons per site, 
where fi^ = fJ^ /V and \x~\ denotes the integer part of x. For simplicity, we assume 
hi = 1 for all sites so that the boson system is in the Mott insulator (MI) phase with 
1 boson per site. In contrast, the fermions can be in any set of Wannier states, since 
for a vanishing tunneling, the energy of the system does not depend on the distribution 
of the fermions in different lattice sites for the homogeneous optical lattice and will be 
very similar for the different repartitions of fermions in the case of small disorder. 

Assume now that the boson-fermion interaction is turned on positive {U > 0) and 
define a = U/V. The presence of a fermion in site i may expel s = or 1 boson 
depending on the interaction strength. The on-site energy gain in expelling s boson 
from site i is AEi = ys(s — 1) — Us + /j^s. Minimizing AEi versus s, we find that 
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a — /ip +1 that we assume to be 1 in all sites. Within the occupation number 
basis, excitations correspond to having 1 boson in a site with a fermion, instead of 
boson and, therefore, the corresponding excitation energy is ~ U. In the following, we 
assume that the temperature is smaller than U so that the populations of the above 
mentioned excitations can be neglected. 

Summarizing the discussion above, we have assumed that < /ip < V (so that 
hi = 1), U — V < jj}* < U (so that in the lowest energy states, each lattice site is 
populated by either one boson or one fermion but never or 2 particles) , and k^T < U 
(so that only these lowest energy states are significantly populated). The physics of 
our Fermi-Bose mixture can thus be regarded as the one of composite particles made 
of one fermion plus one bosonic hole. Notice that the bosonic hole is created because 
the initially present boson is expelled from the site by the fermion, due to the repulsive 
interactions between bosons and fermions. Within the picture of composites, a lattice 
site is either populated by one composite (i.e. by one fermion plus bosonic hole) or free 
of composite (i.e. populated by one boson). In particular, the vacuum state corresponds 
to the MI phase with 1 boson per site. The annihilation and creation operators of the 
composites are pilj : 

F^ = h\f{P (2) 

Fl = Vflh (3) 

where V is the projector onto the sub-Hilbert space of the composites. It is 
straightforward to show that Fi and are fermionic operators. 

2.3. Second- order perturbation: effective Hamiltonian 

The sub-Hilbert space of the composites is [N\/{N^)\{N — A^£)!]-dimensional, a number 
that corresponds to the number of possibilities to distribute the fermions into the 
N lattice sites. We assume now that the tunneling rates J-^ and Jp are small but finite. 
For weak enough disorder, one can assume site independent tunneling rates for bosons 
and fermions ^1]. Using second order projection perturbation theory we derive 
an effective Hamiltonian for the fermionic composites by means of a unitary transform 
applied to the total Hamiltonian (P) . The general expression of the effective Hamiltonian 
can be written as: 



- ^{out\VHMQ 



(owt|ifgg|m) = {out\Ho\in) + {out\VHintP\in) (4) 

rr \ + rr \ QH.ntV\in) . 

where Q = 1 — V is the complement of the projection operator V. As second order 
theory can only connect states that differ on one set of two adjacent sites, -ffgff can only 
contain nearest-neighbor hopping and interactions as well as on-site energies /Zj pTSlj : 

^eff = E [-d^^F^F^ + ^■^] + E K^,M,M, + T^.M, (5) 
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where Fi and F- are defined in Eqs. (jSHSI) and Mi = F^Fi is the composite number 
operator in lattice site i. Ahhough the general form of Hamiltonian (jSj) is universal 
for all types of composites, the explicit calculation of the coefficients dij, Kij and /Zj 
depends on the concrete type of composites that can be created for different values and 
sign of alpha j2Sl- The nearest neighbor hopping for the composites is described by 
—dij and the nearest neighbor composite-composite interaction is given by Kij, which 
may be repulsive (> 0) or attractive (< 0). This effective model is equivalent to that of 
spinless interacting fermions [S]. 

In the situation under consideration, each site contains either one boson or one 
fermion. Therefore, a fermion jump from site i to site j can only occur if the boson 
that was initially in site j jumps back to site i into the hole the fermion leaves behind. 
Therefore, the number operator for fermions and bosons are related to the number 
operator of composites, i.e. = = 1 — n^. From Eq. (@]), we find 



d. 



K - = ^ 




(6) 



(7) 



b 



V ^ 



(id) 

1 



a + aI 



(8) 



with A 



f,b 



^f,b 



^f,b 



Hj' measures the inhomogeneity of the lattice sites. Here, {i,j) 
represents all sites j adjacent to site i. The coupling parameters in Hamiltonian (0) 
depend on all parameters (Jj^, Jf, U, V, /if, jj,}) of the FBH model ((H) in a complex 
fashion. For example, the hopping amplitudes dij depend on disorder but always remain 
positive. In contrast, the interaction term Kij can be positive or negative. All terms 

' b f 

dij, Kij, /ij are random due to the inhomogeneous on-site energies /if and /i- of the 
bosons and the fermions respectively. It is thus clear that the qualitative character of 
interactions may be controlled by the inhomogeneities ^2)- 

For example, consider the case where disorder only applies to the bosons (/i- = est, 
uniform) |32]. The corresponding coupling parameters are plotted in Fig. ^ As 
mentioned above, the hopping amplitudes dij are always positive, although may vary 
quite significantly with disorder, especially when a|^- ~ a. As shown in Fig. Q for 
a > 1, Kij < and we deal with attractive (although random) interactions. For a < 1, 
Kij > and the interactions between composites are repulsive. For a < 1, but close to 



Kij might take positive or negative values for A|j- small or a|^- 



a. 
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Figure 1. (color online) Tunneling, dij, and nearest neighbor interactions Kij between 
the fermion composites as a function of the disorder of bosons aJ^ for various boson- 

r 

fermion interactions a. The disorder for fermions is assumed to vanish here {fj.\ — 0). 



In the following two sections, we investigate the ground state properties of the 
Fermi-Bose mixture (or equivalently that of the fermion composites) in the presence of 
diagonal disorder. We distinguish two limiting cases. The first corresponds to small 
disorder while the second corresponds to A^-''^ ^ a that interestingly maps to a spin 
glass problem. 

3. Quantum phase diagrams of Fermi-Bose mixtures in weakly disordered 
optical lattices 

The first limiting case corresponding to weak disorder (A}y ^ is investigated 
in this section. For the sake of simplicity, we use J\) = = J ■ In this limit, the 
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contributions of the inhomogeneities A^y to the tunnehng dij and to the interaction Kij 
terms in the effective composite Hamiltonian (jSj) can be neglected {dij ~ d and ~ K) 
and we keep only the leading disorder contribution in Jl^ [first term in Eq. (jH))]. Note, 
that the latter contribution is particularly relevant in ID and 2D leading to Anderson 
localization of single particles [33^. We will first describe the non-disordered case and 
then we will discuss the effect of weak disorder. 

The analysis reported below is supported by mean-field numerical calculations (for 
details about the numerical method, see in Ref. jlEl)- We consider a 2D optical lattice 
with = 100 sites with A^^ = 60 bosons, A^f = 40 fermions and J^/V = J^/V = 0.02 
to compute the ground state of the system in the presence of a very shallow harmonic 
trapping potential = u^^^ x l{iY, where l{i) is the distance from site i to the center 
in cell size units) with eventually different amplitudes for bosons and fermions. This 
simulates optical or magnetic trapping which turns out to be hardly avoidable in current 
experiments on ultracold atoms. In addition, it breaks the equivalence of all lattice sites 
and makes more obvious the different phases that one can achieve (see below). In the 
numerics, we have used u = and u = 5 x 10"''. 

In the absence of interactions between the bosons and the fermions {a = 0), the 
bosons are well inside the MI phase with n = 1 boson per site in the center of the trap 
[3 EH] • Besides, the non- interacting fermions are delocalized and due to the very small 

r 

value of u they do not feel significantly the confining trap as shown in Fig. |2ta). 

3.1. Quantum phases in non- disordered optical lattices 

In the absence of disorder, the physics of the Fermi-Bose mixture is mainly determined 
by the ratio K/ d and the sign of K where K and d are site independent parameters. Once 
the fermion composites are created {a > pP, see section 1221), "^^ have K/d = —2{a — 1) 
as a result of Eqs. (jHllZj) with J-^ = = J and Ajj = 0. We now discuss the quantum 
phases that are accessible depending on the control parameter a ^24^. 

For /i^ < a < 1, the interactions between composite fermions are repulsive and of 
the same order of magnitude as the tunneling {K ~ d). Therefore, the system enters 
an interacting Fermi liquid quantum phase [see Fig. Efb)]. The fermions are delocalized 
over the entire lattice but populate preferably the center of the confining trap. Small 
repulsive composite interactions tend to flatten the density profile compared to that of 
non-interacting composites (see Fig. |21and text below). 

For a ~ 1, although the interactions between the bosons and the fermions are 
large, the interactions between the fermion composites vanish and the system shows 
up properties of an ideal Fermi gas [see Fig. Efc)]. Again, the fermions are delocalized 
and their distribution follows the harmonic confinement. 
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Figure 2. Local on-site probability of finding one fermion and zero boson at each 
lattice site for the N = 100 lattice sites in a Fermi-Bose mixture with N-^ = 60, 
= 40, Ji)/V = Jf/V — 0.02 and in the presence of harmonic traps for bosons and 
fermions characterized by cj^ = 10^^ and uj^ = 5 x 10^^, respectively. The interaction 
between fermions and bosons is (a) a = [independent bosonic MI and Fermi gas], (b) 
a = 0.5 [Fermi liquid], (c) a = 1 [ideal Fermi gas] and (d) a = 10 [fermionic insulator 
domain] . 

Growing further the repulsive interactions between bosons and fermions, the 
interactions between the fermion composites become attractive. For 1 < a < 2, one 
expects the system to be a weakly interacting superfluid, whereas for a > 2 a fermionic 
insulator domain phase is predicted [see Fig.|21^d)]. In this case, the fermions are pinned 
in the lattice sites. They tend to merge because of site-to-site attractive interactions 
and populate the center of the trap. 
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Notice that contrary to the bare fermions, the composite fermions are significantly 
affected by the harmonic trapping potential. This is because the coupling parameters of 
the composite effective Hamiltonian ^ are much smaller than the coupling parameters 
of the bare fermion-boson Hubbard Hamiltonian (^. Therefore the harmonic potential 
is able to compete with tunneling and interactions for the composites |2(ij . 

3.2. Quantum phases in disordered optical lattices 

We now assume that small on-site inhomogeneities are present and we investigate the 
effect of disorder on the quantum phase diagram of the system depending on the 
parameters of the effective Hamiltonian (jSJ. Fig. El shows a schematic representation of 
expected disordered phases of the fermionic composites for small disorder. 



Domain Fermi Mott Insulator 

insulator .Dirty" Glass jj; (Checker-board 
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>1 and <0 







»1 and > 



^d 



Figure 3. (color online) Schematic phase diagram of fermionic composites for small 
disorder (aJ^- <g; 1, a) as a function of the ratio between nearest neighbor interactions 
and tunneling of the composites [K/d). 



For <^ d, the system is in the Fermi glass phase, i.e. Anderson localized (and 
many-body corrected) single particle states are occupied according to the Fermi-Dirac 
statistics [32]. The effect of disorder is to localize the fermions preferably into the 
deepest sites. 

For repulsive and large composite interactions (K > and ^ rf), the ground 
state is a Fermi MI phase and the composite fermions are pinned at large filling factors 
preferably into the deepest wells. At half filhng factor (p^ = 1/2), one expects the 
ground state to form a checker-hoard., i.e. the lattice sites are alternatively empty and 
populated by one composite. 
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For large attractive composite interactions {K < and li^l ^ d), the fermions form 
a domain insulator which average position results from the competition between the 
random and the confining potential. 

Finally, for intermediate values of K/d, with K > 0, delocalized metallic phases 
with enhanced persistent currents are possible jSHj- Similarly, for attractive interactions 
{K < 0) and \K\ < d one expects a competition between pairing of fermions and disorder 
i.e. a dirty superfluid phase. 

Further information is provided by numerical computations that we present now. 
We consider on-site random inhomogeneities for the bosons We start from a non- 
disordered phase [A(t = 0) = 0] and we slowly increase the standard deviation of the 

disorder '\/(G^F)^)"^-Tc^F))^ = '^i^) from to its final value A. 

We first study the transition from a (composite) Fermi gas in the absence of disorder 
[see Fig. IH^b)] to a (composite) Fermi glass [see Fig. life)]. Initially [A(t) = 0], the 
composite fermions are delocalized although confined near the center of the effective 
harmonic potential [{u^ — Uf,) x /(z)^]. The local populations fluctuate around {mi) ~ 0.4 
with a standard deviation a/ {{mi — {mi)Y) ~ 0.43. Increasing the amplitude of 
disorder, the fluctuations of decreases as shown in Fig. Ufa). This indicates that 
the composites localize more and more in the lattice sites to form a Fermi glass. For 
A = 5 X 10~^, the composite fermions are pinned in random sites as shown in Fig. EJj^c). 
As expected, the composite fermions populate the sites with minimal jif . 

We now turn to the transition from a Fermi insulator domain phase [see Fig. Efb)] 
to a disordered insulating phase while slowly increasing the amplitude of the disorder. 
For the Fermi insulator domain, the composite fermions are pinned near the center of the 
harmonic trap and surrounded by a ring of delocalized fermions which results in finite 
fluctuations on the fermion occupation number {^J {{mi — {mi)Y) ~ 0.35). As shown in 
Fig. Efa), while ramping up the amplitude of disorder, the fluctuations decrease down 
to ^J {{mi — {mi)Y) ^ 0.13 for A > 10^'^ showing that the composite fermions are 
pinned in different lattice sites. This can be seen in the plot the site population of the 
composite fermions presented in Fig. Efc). Contrary to what happens for the transition 
from the Fermi gas to Fermi glass, the composites mostly populate the central part of 
the confining potential. This is because (i) the attractive interaction between composites 
is of the order of ~ —1.4 x 10"'^ and competes with disorder (A = 3 x 10~^) and (ii) 
because tunneling is small (rf ~ 8 x 10~^) so the fermions can hardly move during the 
ramp of disorder. 

4. The spin glass limit 

4.1. From composite Hamiltonians to spin glasses models 

The second limit corresponds to a case where the interactions between fermions and 
bosons are of the order of, but slightly smaller than the interactions between the bosons 



Strongly correlated Fermi-Bose mixtures in disordered optical lattices 



12 




0-0 1.0x10"^ 2.0x10^ 3.0x10"^ 4.0x10"^ 5.0x10"^ 

A 




Figure 4. Dynamical crossover from the Fermi gas to the Fermi glass phases. The 
parameters are the same as in Fig. |2Ic) . (a) Variance of the number of fermions per 
lattice site as a function of the amplitude of the disorder A. (b) Probability of having 
one composite (one fermion and zero boson) at each lattice site for the M sites in 
the absence of disorder and (c) after ramping up adiabatically diagonal disorder with 
amplitude A = 5 x 10^"^. 
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Figure 5. Dynamical crossover from tlie fermionic domain insulator to a disordered 
insulating phase. The parameters correspond to Fig. Efd). (a) Variance of the 
number of fermions per lattice site as a function of the amplitude of the disorder, 
(b) Probability of having one composite (one fermion and zero boson) in each lattice 
site in the absence of disorder and (c) after ramping up adiabatically the disorder with 
amphtude A = 3 x 10"''. 
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(a ~ 1 and a < 1). As shown in Fig. ^ the effective interaction term Ki j has a zero 
at large inhomogeneities and varies strongly with |A|^j| reaching both positive 

to negative values. Such a situation is accessible for ultracold atomic systems using a 
superlattice with a spatial period twice as large as the lattice spacing plus a random 
potential, both acting on the bosons. The interaction with this superlattice results in an 
alternatively positive and negative additional on-site energy of the bosons, and whose 
amplitude is controlled by the intensity of the superlattice. In particular, one can set 
|Aj jl ~ a. An additional weak random potential introduces disorder. 

Due to the random on-site effective energy /Jj, the effective tunneling becomes non- 
resonant and can be neglected in first approximation while Kij is random with a given 
average (eventually zero) and strong fluctuations from positive to negative values. The 
effective Hamiltonian Q then reduces to 

where we have introduced Si = 2Mj — 1 = ±1. Interpreting Si as classical Ising spins |H7j . 
this Hamiltonian is equivalent to the well known Edwards- Anderson model [^Hl- This 
describes spin glasses, i.e., a Ising model with random positive (anti- ferromagnetic) or 
negative (ferromagnetic) exchange terms Kij. Our system however differs from the usual 
Edwards- Anderson spin glass model as (i) it has a random magnetic term /Zj and (ii) the 
average magnetization per site m = 2N^/N — 1 is fixed by the total number of fermions 
in the lattice. It however shares basic characteristics with spin glasses as being a spin 
Hamiltonian with random spin exchange terms Kij. In particular, this provides bond 
frustration, which is essential for the appearance of the spin glass phase and turns out 
to introduce severe difficulties for analytical and numerical analyses. We thus think that 
experiments with ultracold atoms can provide a useful quantum simulator to address 
challenging questions related to spin glasses such as the nature of the ordering of its 
ground- and possibly metastable states jH EHl ED] , broken symmetry and dynamics of 
spin glasses jHl El] ■ 

In the following two sections, we outline some general properties of spin glasses 
and then we apply the replica method under the constraint of a fixed magnetization and 
argue that this preserves the occurrence of a symmetry breaking characteristics of spin 
glasses in the Mezard-Parisi theory |4:. 

4-2. Generalities on spin glasses 

Consider a spin glass at finite temperature with a random exchange term K^j with 
average K and variance AK. The magnetization is characterized by two order 
parameters: (i) m = (sj), the average magnetization per site and (ii) ^ea = (si)'^, 
the Edwards-Anderson parameter, where ~ denotes the average over disorder while (■) 
represents the thermodynamics average. It is clear that m ^ signals a long-range 
magnetic order while ^ea 7^ signals a local magnetization that may vary from site to 
site and from one configuration of quenched disorder to another. 
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Figure 6. Schematic ptiase diagram of a spinglass. 



Earlier experimental studies have identified three magnetic phases as schematically 
represented in Fig. IHl : (i) At high temperature and small average spin exchange K, 
one finds a paramagnetic phase characterized by m = = 0. (ii) For > and 
large, one has a ferromagnet with m 7^ and ^ea 7^ 0. (iii) For weak K and small 
temperatures, a spin glass phase appears with m = but ^ea 7^ 0. This signals that the 
local magnetization is frozen but that disorder prevents a long-range magnetic order. 

As pointed out before, the physics of spin glasses still opens challenging questions. 
In particular, there are competing theories that predict different natures of the magnetic 
order. 

The 'droplet' picture is a phenomenological theory based on numerical and scaling 
arguments. It predicts that there are two ground states related by spin-fiip symmetry 
and that low-lying excitations are domains with fractal boundaries (the droplets) with 
all spins inversed compared to the groundstate. This theory is supported by numerical 
computations and is believed to be valid in short-range spin glass models such as the 
Edwards- Anderson [5H] . 

The Mezard-Parisi picture predicts a large number of low-energy states with very 
similar energies. This leads in particular to disorder-induced quantum frustration. This 
theory is a mean-field, based on the replica method which is a special trick introduced 
to compute the non-trivial average over disorder of the free energy functional. The 
Mezard-Parisi theory is formulated in long-range spin glasses such as the Sherrington- 
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Kirkpartick model 

1 V 1 . 

^S-K = 4 ^ij^^^J + 2 

that differs from the Edwards-Anderson model (jH)) by the long-range spin exchange 
[{ij) denotes here the sum over all pairs of lattice sites, either neighbors or not]. The 
applicability of the Mezard-Parisi picture to short-range spin glasses is still controversial. 



4-.3. The replica- symmetric solution for fixed magnetization 



As the disorder is quenched, one must average over disorder the free energy density, 
/ = — In Z/k-Q using the replica trick. In the following we assume that the random 
variables Kij and /Ij are Gaussian distributed with average K and h respectively and 
standard deviation AK and Ah respectively. We form n identical copies of the system 
(the replicas) and the average is calculated for an integer n and a finite number of spins 
A^. Then, using the well-known formula In x 
the analytic continuation of for n 
N ^ oo. Explicitly, is given by: 



\imn^o(x^ — l)/n, InZ is obtained from 
0. Finally, we take the thermodynamic limit 



= ^ exp -H[s 



n 



'111 



=±1 



where H[sf,n\ is the sum of n independent and identical spin Hamiltonians, averaged 
over the disorder, with Greek indices now numbering the n replicas. Computing 
the average over disorder leads to couphng between spin-spin-interactions of different 
replicas. 

After some analytics that are detailed in one finally gets 



/SK _ (K/kj^TY^ 
k^T 4 



{Ah/k^TY^ 



(12) 



+ / dz 



In 



2 cosh 



^K^q + Ah^ - h 



and 



m 



dz 



dz 



tanh 



tanh 




(13) 



(14) 



These are characteristic values of spinglasses that may be measured in experimental 
realizations of the proposed systems. 

Finally, a study of stability of the replica method as detailed in [211 shows that 
the magnetization constraint specific to our model would not change the occurrence of 
replica symmetry breaking, provided the Mezard-Parisi approach is valid in finite range 
spin glass models. 
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5. Conclusion 

In this paper, we have reviewed our recent theoretical works on Fermi-Bose mixtures 
in disordered optical lattices. In the strongly correlated regime and under constraints 
that we have discussed in detail, the physics of the mixture can be mapped into that of 
single-species Fermi composites. This is governed by a Fermi-Hubbard like Hamiltonian 
with parameters that can be controlled with accuracy in state-of-the-art experiments 
on ultracold atoms. 

We have shown that the presence of disorder (created by random on-site energies) 
introduces further control possibilities and induces an extraordinary rich quantum phase 
diagram. For the sake of conciseness, we have restricted our discussion to a particular 
regime that proves to be a case study (more details may be found in Ref. [2n])- For 
weak disorder, we have discussed the phase diagram which corresponds to the formation 
of Fermi glass. Domain insulator, dirty superfiuids and metallic phases. Numerical 
calculations support our discussion. 

For larger amplitudes of disorder, we have shown that the Hamiltonian reduces to 
that of a spinglass, i.e. a spin system with random exchange terms. In our system, 
the fictitious spins are coded by the presence or the absence of particles in each lattice 
site. The physics of spinglasses is a challenging problem in statistical physics which 
is still unsolved. In particular, two theories are competing: the droplet picture and 
the Mezard-Parisi picture. These two theories lead to different predictions and even 
the nature of ordering in spinglasses is not known. On the one hand, the Mezard- 
Parisi picture, which is assumed to be valid for long-range spin exchange, predicts 
the existence of a huge number of quantum states with very similar energies that all 
contribute to the low-temperature physics of spinglasses. On the other hand, the droplet 
picture, which is believed to be valid for short spin exchange, assumes the existence of 
only two pure states connected by spin flip symmetry and excitations are domains of 
constant magnetization with fractal boundaries. As possible experimental realizations 
of spinglass systems with controllable parameters, mixtures of ultracold fermions and 
bosons may serve as quantum simulators to solve the controversy and shed new light on 
this extraordinary rich physics. 
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